#ifndef SIMULATION_HPP
#define SIMULATION_HPP
#include "mediator.hpp"
namespace QuantumWalk{
  template<typename T>
  class simulation{
    
  private:
    int nstep_;
    mediator<T>* p_med_; 
  public:
    simulation();
    simulation(int nstep,Imediator* p_med):
      nstep_(nstep){
      p_med_=static_cast<mediator<T>* >(p_med);
    };
    //vector<double> p_sites(){};
    void step(){};
    void run(){};
    //vector<double> avg_p_sites(){};
    //double mixing_time(vector<double> avg_p_full,double threshold){};
    //double sampling_time(vector<int> sites,vector<double> avg_p,
    //			 double threshold){};
  };
};
#endif
